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I. INTRODUCTION 



Polymeric fluids exhibit a rich and complex set of phenomena associated both with system-specific and 
global properties of the polymer molecules. Chemical details become increasingly less relevant for global 
polymer properties as the degree of polymerization N increases J^i^i^i^ Thus, for TV oo, one can use 
coarse-grained models in which only the most fundamental aspects of the polymer structure are taken into 
account. The behavior of polymer solutions depends in general on temperature. For T large enough, the 
most relevant feature is the local repulsion. In this regime, usually called good-solvent regime, the radius of 
gyration Rg, as well as any other quantity that is related to the global size of the polymer, scales as N'^, 
where is a universal exponent; v w 0.5876 (Refi^). As T is lowered, one reaches the 6 temperature, Tg, 
below which polymers are compact (Rg ~ N^^^) and phase separation occurs.^' At the 6 point polymers 
behave approximately as Gaussian coils. The crossover from good-solvent to behavior is well understood. 
For N ^ oo any global quantity O behaves as)^ 

0{T,N,c) = aiOG{N,c)fo[a2{T -Te)N^^^\nN)-^/^\a3cRl{T,N)]. (1.1) 

Here c is the polymer number density, Rg{T, N) is the zero-density radius of gyration, and Og{N, c) is the 
expression of O for ideal chains. The function fo{x,y) is universal, all chemical details being included in 
the constants a^. 

Eq. (iri|) is strictly valid only for iV oo, T -> Tg, at fixed a2{T -Tg)N^/'^{\-RN)-'^l'^^ . For finite values of 
N one should also take into account the corrections to Eq. p.ip that decay very slowly, as inverse powers of 
InA'^. In Ref.— we computed the crossover curve for the interpenetration ratio 5* and found that logarithmic 
corrections are only relevant very close to Tg. Outside a tricritical region around Tg, Eq. (jl.ip provides a 
reasonably accurate description of the crossover. As emphasized in Rcf.^*', in order to compute the crossover 
functions defined in Eq. (jl.ip one can use the continuum two- parameter model (TPM)fii Indeed, if we 
identify {T — Tg)N^/'^ {In N)~'^/^^ with the Zimm-Stockmayer-Fixmani^ variable z (with a model-dependent 
proportionality factor), then the crossover function for O corresponds exactly to its TPM expression. 

It is interesting to note that the TPM is also of interest to describe the corrections to scaling in some 
polymeric systems. Indeed, as discussed in Ref J^, the TPM describes the approach to the scaling limit when 
Vm/l^ ^ 1, where I is the persistence length and Vm is the volume occupied by a polymer blob of length I. 

In this paper we wish to compute the crossover functions for several quantities whose behavior in the 
good-solvent regime has been considered in RefsJ^ii^. We compute numerically the TPM predictions for the 
second, third, and fourth virial coefficient, for the swelling factors, and the density corrections to the radius 
of gyration, the end-to-end distance, and the hydrodynamic radius. This allows us to give exact predictions 
for the thermodynamic behavior and for the polymer size in the whole dilute regime $p ^ 1, where $p is 
the polymer volume fraction. 

The paper is organized as follows. In Sec. |TT]we define the TPM as the scaling limit of the lattice Domb- 
Joyce model. This is a rigorous well-defined definition that does not rely on perturbative field theory. In 
Sec, mil we define the quantities that are considered in the paper and report some results and properties that 
are useful in the following section. In Sec. IIVI we give the results of our work. We first report the analysis 
of the Monte Carlo results and then determine the TPM functions associated with the different quantities. 
In Sec. IIV CI and IIVDI we use these results to predict the osmotic pressure and the polymer size in the 
dilute and in the semidilute regime. Finally, in Sec. IIV El we compare our predictions with the available 
renormalization-group results. Some conclusions are presented in Sec. IVl 



II. THE DOMB-JOYCE MODEL 



In order to compute the TPM crossover functions, we consider the three-dimensional lattice Domb- Joyce 
(DJ) modelji^ We consider a cubic lattice and model a polymer of length iV as a random walk {tq, ri, . . . , ttv} 
with \ro> — fa+il = 1 on a cubic lattice. To each walk we associate a Boltzmann factor 

^-pH^^-n^a^ J2 5r^,,^, (2.1) 

Q<a<(3<N 

with w > 0. The factor a counts how many self-intersections are present in the walk. This model is 
similar to the standard self-avoiding walk (SAW) model, in which polymers are modelled by random walks 
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in which self-intersections are forbidden. The SAW model is obtained for w = +00. For finite positive w 
self-intersections are possible although energetically penalized. For any positive w, this model has the same 
scaling limit of the SAW modeli^ and thus allows us to compute the universal scaling functions that are 
relevant for polymer solutions. 

The DJ model can be efficiently simulated by using the pivot algorithmj ^^'^^i^^ i?° For the SAW an efficient 
implementation is discussed in Ref.^^. The extension to the DJ model is straightforward, the changes in 
energy being taken into account by means of a Metropolis test. Such a step should be included carefully 
in order not to loose the good scaling behavior of the CPU time for attempted move. We use here the 
implementation discussed in Ref.^^. 

The TPM results can be derived from simulations of the DJ model. Indeed, the continuum results are 
obtained^^ by taking the limit ^ 0, A'' ^ 00 at fixed product wN^^"^ (we call it x). The variable x 
interpolates between the ideal-chain limit {x — 0) and the good-solvent limit {x = 00). Indeed, for w = 
the DJ model is simply the random-walk model, while for any w ^ and N 00 one always obtains the 
good-solvent scaling behavior. The variable x is directly related to the variable z that is usually used in 
the TPM context'!^ indeed, z — ax. The normalization factor a can be fixed by considering the small-z 
behavior of the interpenetration ratio conventionally one takes 'i' = z + O(z^). In the DJ model^i 
\E' = (3/27r)^/^i(;Afi/2 for small wA^^/^, so that we can identify 



III. DEFINITIONS 



We consider the osmotic pressure n(T, A^, c) or, equivalently, the adimensional compressibility factor 

where c is the polymer number density, p the weight concentration, M the molar mass of the polymer, T 
the absolute temperature, ks and R the Boltzmann and the ideal-gas constants. In the dilute limit, Z can 
be expanded in powers of the concentration as 

z = i + Y,McRir-\ (3.2) 

n>2 

where Rg the zero-density radius of gyration. The coefficients An depend on T, N, and on chemical details. 
However, in the good-solvent regime, they converge to universal constants ^* as A^ — > 00. Moreover, the 
renormalization group predicts that corrections should always scale as N^^, where A is a universal exponent 
whose best estimate is^ A = 0.515 ± 0.007t^;^JEJ. Therefore, for large A^ we expect 

A„(r, N) - Al + Ai,n{T)N~^ + ■■■ (3.3) 

While the constants A*^ are universal, the coefficients ^i,„(T) are system specific and temperature dependent. 
However, the ratios 



_ A,AT) 

are also universal. Precise estimates of A„ for n < 4 and of 63 have been obtained in RefiM: 



(3.4) 



^2 = 5.500 ±0.003, (3.5) 

^3 = 9.80 ±0.02, (3.6) 

Al = -9.0 ±0.5, (3.7) 

63 = 4.75 ±0.30. (3.8) 

Instead of A2 it is customary to define the interpenetration ratic^^ 5* = 2{Att)^'^/'^ A2, whose large- A^ value 

in the good-solvent regime is^'^ = 0.24693 ± 0.00013 (Ref.^"*), = 0.24685 ± 0.00011 (Ref,^). 
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In the TPM, the coefRcients An become functions of z, An{z), such that A„(z = 0) = (z = is the 
ideal-chain case) and A„(z — oo) = A* . Since z ^ N^^^, Eq. (13 3p imphes 

A„(z) = A; + a„z-2A (3.9) 

for large z, with a„/a2 = 6™.-^ 

The small-z behavior of the TPM functions can be determined by using perturbation theoryi^iii We quote 
here the result for ^2(2) and ^3(2)1^1^1^ 

A2{z) = i(47r)3/2z[l- 4.779663z+ 25.58964^2 + 0(z3)], (3. 10) 

A^iz) = ^(^\ (208V2-108V3-103)z'VO(z'*)~1100.7z^ (3.11) 
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In the appendix we compute the leading contribution to A^lz) obtaining 

A4(z) = (14075 + 12624\/2 - 18468\/3)7r9/2z4 + ©(z^) w -29883.1zl (3.12) 

^ ' 45045 

Beside the osmotic pressure we consider three different quantities that characterize the polymer size: the 
radius of gyration Rg, the hydrodynamic radius Rh, and the end-to-end distance. In the DJ lattice model 
they are defined as follows: 

Rl = ((ro - VNf) . (3.15) 



We also define the ratios 



A,H ^ 1^, (3.16) 



Rl Rh 
where a hat indicates a zero-density quantity, and consider the density expansions 

|| = l + S^,g{cRl)+S2,g{cRlf + ■■■ 



= 1 + Si.AcRl) + S2A(=Rl) 



^ = 1 + Si.h{cRI) + S2M^RI? + ■■■ (3.17) 
Rh 

The ratios (j3.16p and the density coefficient Sn,=if: are system-dependent quantities. However, as — > cxd in 
the good-solvent limit, they approach universal quantities, which will be labelled as A*g^, ^gff' ^^'^ ^n#- 
The limiting values of the ratios and of the density coefficients for n = 1,2 have been determined in Reff^^. 

In the TPM all previous quantities are functions of z which converge to their good-solvent value for z ^ 00. 
If Q corresponds to Age or to a coefficient Sn for the radius of gyration and the end-to-end distance, we can 
also determine the corrections for z — > cx). Indeed, in this limit, we have 

Q^Q* + aQZ~^^. (3.18) 

The ratio aQ/a2 [02 is defined in Eq. (|3.9p ] is universal; estimates for Age, Si^g, and Si^e are reported in 
Reffi^. For z ^ 0, Sn,i^{z = 0) = 0, while Age and AgH converge to the ideal-chain (random-walk) values 

Age{z^O) = ^ Agnize 0) = ^. (3.19) 
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The leading corrections for z — > to all these quantities are reported in the Appendix. 
Finally, we consider the swelling factors for the zero-density radii: 



= ^Ne^aliz), (3.20) 

Rl = Nfaliz), (3.21) 

1 /o \ 1/2 
1 / OTT ^ ' 



i?H - g I^tJ ^N^(^h{z) . (3.22) 

The swelling factors are normalized so that a^{z = 0) = 1. In the DJ model the metrical factor £ is equal 
to the lattice spacing. For z ^ oo they behave as 

a# = ao,#z^''-\l + ai,#z-2A + . . .) (3,23) 
For the hydrodynamic radius one should additionally consider corrections proportionalSi to z'*""^. 



IV. CROSSOVER FUNCTIONS 
A. Monte Carlo results 



The main purpose of the present paper is the determination of the crossover functions for the quantities 
defined in Sec. lIIII We consider five different values of z, which we denote by zi, . . ., Z5, which belong to the 
crossover region between ideal and good-solvent behavior. Explicitly we use zi = 0.056215, Z2 = 0.148726, 
Z3 = 0.32165, Z4 — 0.728877, and Z5 = 2.50828. They were chosen so that A2(z„) « n (remember that A2{z) 
varies between and 5.50). In order to compute the TPM value for each z^, we perform several simulations 

at values (w^ , iVy) such that WijNl-"^ — {2-k /i)'^/'^ Zi, choosing Nij between 100 and 8000. Then, we fit each 
universal quantity Q with the theoretically expected behavior j^i^i^ 

Q{w,,,N,^) = Q*(zO + iV^'%(z,) + iV^icQ(z,). (4.1) 

The TPM result corresponds to the leading term Q*{zi). In order to detect additional scaling corrections 
that are not taken into account by the fit ansatz (|4.ip . we have repeated the fit several times, each time 
including only data satisfying N > N^i^. 

We illustrate the procedure by considering A2. In Fig. [1] we plot A2{w,N) vs N~^/^; for each z we also 
plot the function Q* + N'^/'^bq obtained in the fit of aU data (iVmin — 100) to Eq. (|4.ip . The data points 
follow the expected behavior quite precisely, with very small N^^ corrections. Note that an extrapolation 
is always needed except for very small values of z. Estimates of A2{z) are reported in Table [J for different 
values of A^min- No systematic deviations are observed. We take the results corresponding to Nn^^ = 500 as 
our final results. They are reported in Table HTl We have applied the same analysis to A3, A^, Age, AgH, 
S7i,g, Sn,e, Sn,H [n — 1, 2). The results corresponding to iVmin = 500 are reported in Table [Til 

Finally, we determine the swelling factors. We consider 6R^/N, R^/N, and hnVN/Rn [kn = 
2-7/2(3^-^1/2 -J, The corresponding quantities are reported in Fig. O Their behavior with N at fixed z 
is perfectly consistent with Eq. (|4.ip . Therefore, we have performed the same fits as before. The results are 
reported in Table HIl 



B. Interpolation formulas 

We now use the results of Sec. lTV A[ the good-solvent results of Refs.— i^^, and the small-z results mentioned 
in Sec, mil and in the Appendix, to obtain interpolation formulas that are valid for all values of z. We discuss 
in detail the virial coefficients; all other quantities are analyzed analogously. 

For the second virial coefficient we wish to find an interpolation that satisfies the following properties: (i) for 
z -> cx) it must satisfy A2{z) ^ A'^ = 5.500 [Eq. ^^]; (u) for z ^ it must behave as 47r3/2z(l -4.779663z) 
[Eq. p.lOp ]: (iii) the interpolating curve should assume the values determined numerically and reported in 
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Table HIl Since the results for ^2(2) indicate that this function is monotonic, we take an interpolating 
function of the form 

A2{z) = 47r=^/2z(l + diz + d2Z^ + dsz^ + d^z'^)-'^''^ , (4.2) 

where di are constants to be determined. The constant di can be fixed to obtain the expansion (|3.10p to 
order z^: we obtain di = 19.1187. The constant d^ can be fixed by requiring ^2(2 — 00) — A2, where A2 is 
given in Eq. (|3.5p : this gives d^ = 268.96. Then, we fit 



1-diz- d^z^ = + da^ • (4.3) 



Using the five data reported in Table HIl we obtain d2 — 126.783 and ds — 331.99. The interpolation formula 
is reported in Table Hill For z 00, Eq. (|4.2p gives 

A2{z) = 5.500 - 1.6972/z + 0{z-^) . (4.4) 

This expression is compatible with Eq. p.9p . taking into account that^** 2 A w 1.03. It allows us to estimate 
a2: 02 ~ —1.7. Of course, this is a very rough estimate. A careful determination would require A2{z) for 
much larger values of z and a careful analysis of the corrections to the behavior (|3.9p . Expression (|4.4p agrees 
with the field-theoretical result reported in Ref.^^, which predicts 02 ~ 5.50 x (—0.30) w —1.65. 

We now compare the interpolation formula (|4.2p with similar expressions that appear in the literature. 
We consider the expression reported in Ref.^ (Sec. 15.5.2): 

1,^ ,3/2 0.1825(1 + 2.15z + 0.82i2)-0-236 
- 2(^") ' (l + 1.32z + 0.378i^)o^^a4 ■ (4-5) 

where22, z — 0.182z. Equation (|4.5p has been obtained by using a sophisticated form of rcnormalized one- loop 
perturbation theory. By means of an extensive Monte Carlo simulation Rcf.^^ obtained 

A2{z) = i(47r)3/2z(i + 14.339Z + 60.30^^ + 66.37^)-^/^. (4.6) 



Finally, we quote the field-theoretical expression of Ref.'^'' for e = 1: 

21r; 



A2(z) = 1(4^)3/2 + l (Aln2+^-)(^) 

^ ' 2^ ' I 8(1 -fry) 64 V 6^1 + '// 



4(1 + ?/) 



(4.7) 



with rj = 256z / 53 (this relation is obtained by matching the small-z behavior) . Our result (|4.2p is essentially 
identical to Eq. (|4.6p : differences are less than 0.3%. It is also in very good agreement with the field- 
theoretical result (|4.5p . differences being less than 1.5%. Eq. (|4.7p is worse: the difference is of order 8% for 
z = 0.1 and increases to 12% for large values of z. 

Let us now discuss ^3. We will use an interpolation formula analogous to (|4.2p . setting 

Aj,{z) = 1100.7z3(l -f diz + d2z'^ + dgz^ + diz^)-^/'^. (4.8) 

The prefactor has been fixed by using Eq. (|3.1ip . To determine the coefficients we use a strategy slightly 
different from that discussed for A2, since we only know the leading small-z behavior of A^ and thus we 
cannot fix di by using perturbation theory. Instead, we make use of the results of Ref.^^ to obtain the large-z 
behavior of A3(z). Since^i 63 = 03/02 = 4.75 [see Eq. dH])], Eqs. gU and (HH) give A^iz) « 9.80-8.062/z 
(again we use the approximation 2 A « 1). If we require Eq. I|4.2p to reproduce this expansion, we obtain 
d^ = 594.386 and d^ = 541.906. Finally, we fit the results of Table HIl to determine di and c?2- The resulting 
expression is reported in Table IIIII 

Finally, we consider A^. In this case we do not know 04/02 and thus we use an interpolation formula with 
only three parameters: 

Ai{z) = -29883.1z''(l + diz + d2z'^ + d^z^y^^^. (4.9) 

The prefactor has been fixed by using Eq. (|3.12p . The constant ^3 is fixed by using Eq. (|3.7p . di and d2 by 
fitting the numerical results of Table |TT| The final expression is reported in Table IIIII In Fig. [3] we report 
the crossover functions for A2, A3, A^. They are monotonic and approach the good-solvent value for z > 5. 
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Similar analyses are performed for the density corrections to the radii. In this case, however, the crossover 
functions are not monotonic. For the second density correction, this is evident from the numerical data. 
For instance, S2.g and S2,e are first positive and increasing, in agreement with Eqs. (|A12p . (jA17|) . reach 
a maximum for 0.5 ^ z ^ 1, and then decrease, converging to the good-solvent value which is negative. 
The density coefficient S2,h behaves in the opposite way, but note that, because of its definition, S2.H is 
equivalent in some sense to —5*2,6 and —S2.g- For the first density correction the nonmonotonicity can be 
inferred by using the results of Ref.^'^. If ^i.^ — SI ^{1 + X^z^'^^), we have A#A2/"2 ~ —0.050 for both 
Si^e and Si^g [02 is defined in Eq. ([3^ ]. Using 02 w -1.697 [Eq. (|474l) ]. we obtain: 

51,9 « -0.3152- 0.0049/^2'^, S'l.e « -0.3853 - 0.0059/2^'^. (4.10) 

Thus, the first density correction vanishes for z = 0, then decreases, becomes smaller than the good-solvent 
value, and eventually converges to it from below. This effect is however numerically very small and thus the 
nonmonotonic approach is in practice irrelevant. The nonmonotonic behavior requires interpolation formulas 
slightly different from those used for the virial coefficients. For instance, expressions like (|4.2p cannot change 
sign and thus are unsuitable for S2,#- Our interpolations are reported in Table [TTTl and plotted in Fig. [H 
Note that the interpolations of 82,^ are not very precise in the region 0.5 ^ z < 1, since here the functions 
change their behavior and we do not have enough data points to identify precisely where the functions reach 
their maximum. 

Finally, let us consider the swelling factors. Because of Eq. p.23p . we use an interpolation of the form 

a# = (1 + biz + b2Z^ + b^z'Y^-'-''^/^ (4.11) 

taking^'* v — 0.58758. The coefficient hi is fixed by requiring to reproduce the small- 2: behavior reported 
in the appendix, while 62 and 63 are obtained by interpolating the numerical data. The results are reported 
in Table IIIII and shown in Fig. [S] Note that ae, otg, an behave in a very similar way, differences being 
tiny. For Ug we also show the prediction of Ref,^, ag = {1 + 7.286z + 9Mz'^ f °^''^^^, and that of Ref;^, 
ttg = (1 + 1.32z + 0.37852)° "**, ^here2^ z = 0.182z. The result of Ref.— is perfectly consistent with ours. 
The field-theoretical result [note that in Fig. [5] it can hardly be distinguished from our result for ae(z)] is 
slightly larger (1% at z = 5): differences are mainly related to the different choice of the exponent v. Finally, 
for the end-to-end distance we mention the result of Ref^^i: ofe = (1 + 7.6118z -I- 12.051352^)° °*'^^*^. Again, 
this expression is in perfect agreement with ours. 



C. The osmotic pressure 



Knowledge of the crossover functions for the lowest virial coefficients provides the osmotic pressure in the 
dilute regime in which $p ^ 1, where $p is the polymer packing fraction, 

AnRl inRl Na 

Na the Avogadro number, M the molar mass of the polymer, c and p the number density and the weight 
concentration, respectively. In Fig. [6] we report the compressibility factor Z defined in Eq. (|3.2p for several 
values of z for $p < 1. In this range of concentrations the virial expansion converges quite welU^ and thus 
our interpolations provide accurate estimates of Z as a function of z and ^p. 

In Refii^ it was shown that a resummation of the virial expansion by using the known large-<I>p behavior 
provides a reasonably accurate expression for Z valid in the whole semidilute region. Here we apply the 
same method to the determination of the leading z"^^ correction to the good-solvent value. Since An = 
A* -I- anZ~^^ for large z, we can write 

Z{z, N, c) « Z*{N, c) + z-^^ZiiN, c) . (4.13) 

The functions Z*(N,c) and Zi{N,c) depend on c and N. However, for N 00, the renormalization group 
predicts that they become universal functions of the packing fraction $p, so that 

Z(z, N, c) « Z*{%) + z-^^Zi{%). (4.14) 
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The good-solvent function Z*{^p) is reported in Ref.— . We will now determine the function Zi{^p). For 
this purpose we determine its large-$p behavior. We expect 

$p^oo. (4.15) 

To fix a, we note that 11 is a function of the monomer concentration Cm = cN but not of the degree of 
polymerization for c ^ cxd (and therefore also Cm oo). Hence, cz"^'^$p should be independent of iV, 
once c has been replaced by Cm- This condition implies 

A + 1 , , 

a = ^. (4.16) 

In order to understand the region in the (c, z) plane in which expansion (I4.14p is valid, we must discuss the 
expected scaling behavior in the large-concentration limit for generic values of z. In the TPM the osmotic 
pressure satisfies the general scaling behavior— 

^ cP[ceJ^N^'^,z). (4.17) 



In the limit z — s- the Flory-Huggins theory applies: for large values of c, 11 is proportional to the square of 
the monomer concentration Cm and to the interaction strength w, so that 

^ wcl,t^ - cz(c^3iV3/2). (4.18) 



For large c, the dependence on TV should disappear at fixed Cm and w, so that the relevant scaling variable 
is z/{ci'^N^/'^) = w/{cm£^)- Therefore, we obtain the scaling behavior 

The function fz{x) is finite for a; — > 0, while for x ^ oo^ consistency with (j4.14p implies 

fz{x) « a/x(3''-2)/(3-i)(i + 5^2,-A/(3.-i))^ (4 20) 

The value fziO) can be determined by noting that for z ^ we can write 

Z « zct^N^''^fz{0) = A2cRl = AT:^/^zc{NlQf/^t^. (4.21) 

This implies /z(0) = (27r/3)3/2/2. 

Eq. (|4.19p indicates that, at fixed large z, the compressibility factor shows two different behaviors. If 4>p 
is large but still z ^ c£^N^/^, Z increases following Eq. (|4.14p . If the concentration is further increased, 
the argument of fzix) decreases and eventually Z « zcf^N^/"^ fz{^) = i^/iizag^^p, i.e. Z becomes linear 
in the concentration. It is clear that this second regime cannot be obtained from extrapolations of results 
in the dilute region. We will thus consider only concentrations such that z ^ cl^N^^^, so that we can use 
Eq. (j4T4)l . 

To determine the osmotic pressure for densities in the semidilute regime satisfying z '3> ci'^N^^^ we expand 
■An{z) — A* + An,i/z for z oo. Thus, we obtain for z ^ oo 

Z w 1 + 1.31303$p -I- 0.558533$^ - 0.122455$^ 

+ i (-0.405182$p - 0.459468$^ -I- 0.0507105$^) , (4.22) 

which is consistent with (|4.14p if we approximate 2 A w 1. In Ref.— we determined an interpolation formula 
for the leading term with the correct large-$p behavior. Here we do the same for the correction term: we 
determine an interpolation formula that has the asymptotic behavior (j4.15p for —^ oo with a given by 
Eq. (|4.16p . and agrees with the previous expansion for <I>p — s- 0. A simple expression satisfying these two 
properties is 

- 0.405182$p(l + 2.30016$p + 1.08734$^)0.493 ^ (4 23) 
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Combining this expression with that obtained in Ref.— we obtain for the comprcssibihty factor 

/ 1 + 1.52605*. + 0.795366<J>f \ ^'^^ 

Z = 



1 + 1.52605$p + 0.795366$2 



, 1 + 0.5245$p / 

„M^<i>p(l + 2.30016$^ + 1.08734<i>2)0-493. (4.34) 

It is not possible to determine a priori, for each z, the density range $p < ^p^maxiz) in which Eq. (|4.24p 
apphes. For <i>p = 1 we can compare expression (j4.24p with the virial expansion (|3.2p including the terms 
up to n = 4. The relative difference is less than 5% (1%, 0.1% respectively) for z > 1.7 (4, 12 respectively). 
Thus, for $p = 1, Eq. (|4.24p is substantially correct for 2 > 4 and reasonably predictive for z >2. Accepting 
an error of 5% (it makes little sense to require a smaller error since our interpolation formulas cannot in any 
case be more precise than 5-10%; see the discussion reported in RefJ^), we can set ^p,niax{z = 2) = 1. An 
estimate of $p.max(^) for larger values of z can be obtained by noting that 

$p,max(^) ~ z'''-' ~ z2-53. (4.25) 

Indeed, Eq. (I4.24p is vahd as long as z » ci^N^/'^ ^ cRgag^ ^ <i>pa~^. Hence, since Ug ^ z^"^^ for 
large z, we obtain Eq. (|4.25p . Therefore, with errors at most of 5% we expect the range to extend up to 

%,n.M ~ (z/2)2-53. 

The prediction (|4.24p for $p < 10 is reported in Fig.[7]for several values of z. It is clear that the results for 
z = 2 do not extend beyond $p = 1, since, by increasing $p, Z begins to bend in an unphysical way. No such 
phenomenon is observed for z > 5, which is therefore expected to be the range of z in which (j4.24p applies 
for < 10. This is in agreement with the estimate of 3'p,max(z) given above. Note that scaling corrections 
are quite large in the semidilute regime. For instance, consider z — 10. Since A2/A2 = 5'/^'* = 0.97, in 
the dilute regime the solution is essentially in good-solvent conditions. For $p = 10 we obtain Z = 31.3 to 
be compared with the good-solvent value Zqs — 35.8: the pressure is lower by 14%, a significant deviation 
from the good-solvent value. 

D. Concentration dependence of the polymer size 

The considerations we have presented for the osmotic pressure can be generalized to the radii. Deep in 
the semidilute region, polymers behave like ideal chains and therefore behaves as 



NeUwIc^) = Nffn (^^) , (4.26) 



with /fl(0) ^ 0. As before, we focus on the deviations from the good-solvent regime. If R is the end-to-end 
distance or the radius of gyration, corrections scale as z~^^ (this is not the case for the hydrodynamic radius 
which may show additional corrections proportional toi^i^ z^"^"**). In order to obtain the behavior for large 
c in the good-solvent regime we write 

i?2 = afli?2$^i(l + 6;jz-2A$02)^ (4 27) 

and require this expression to be consistent with Eq. (|4.26p . This allows us to identify ai and a2'- 

2v-l A , , 

ai = -7. 7 "2 - 7. 7- 4.28 

31^—1 31^—1 

Using these expressions, we can now extrapolate our virial results to the whole semidilute region, obtaining 
(the leading terms already appear in Ref J^) 

r)2 1 

-r^ = (l-l-0.801$p + 0.37$2)-o ii5- _0.024$p(l-h0.76$p)-°-55^, (4.29) 

z 

R^ 1 

-J- = (l + 0.655$p-l-0.28$^)~o ii5_ _o.021$p(l + 1.46$p)"°-^5'*. (4.30) 
Rj z 



10 



As discussed before, these expressions are only valid for z <C c^N'^l'^ . We do not present an extrapolation for 
Rh: because of the presence of two different corrections (one proportional to z~'^^ and one proportional to 
^4i/-4^ see Ref.'^^), which make extrapolations of the form (|4.27p incorrect. The results for Rg/Rg are shown 
in Fig.[8]for the same values of z that occur in Fig.[7l Note that for large z, our extrapolation predicts Rg/Rg 
to decrease when z decreases. This is consistent with the nonmonotonic behavior of 6*1,3(2;) we mentioned in 
Sec. IIV Bl and with the numerical results of Reff^^, which suggest a negative scaling correction for S2,g{z): 
S2,g{z) = -0.087 - (0.003 ± 0.005)z"2A gf course, this behavior should eventually change, since R^/Rl 
should converge to 1 as z ^ 0. 

For the radius of gyration corrections are less evident than in the case of Z. For instance, for $p = 10 the 
relative difference between R^/R^ for z = 10 and for z = 00 (good-solvent value) is only of 0.7%, indicating 
that the polymer size is less sensitive to the solvent quality far from the ideal-chain limit. 

E. Comparison with previous renormalization-group results 

In the previous sections we have obtained predictions for the osmotic pressure and the radii in the TPM. 
We wish now to compare these results with those obtained by using field theory. We compare mainly with the 
results reported in Ref,-S=, which have been obtained by using renormalized one-loop perturbative expressions 
with a careful choice of the renormalization constants. In Fig. [H] we report the compressibility factor Z as a 
function of z for $p = 1, i.e. for the largest value of the density at which the virial expansion is supposed to 
work. As we discussed in Ref.-^ this expression should be quite accurate, deviations being at most of 1-2%. 
We also report the field-theoretical result, which is in perfect agreement with our estimate. For larger values 
of we cannot use the virial expansion. Instead, we employ the approximate expression (|4.24p . which is 
valid only for large values of z. In Fig. [TOlwe report Z in the range < $p < 10 for two values of z: z — 9.90 
corresponding toi? = 1 — ^/5'* — 0.03 and z = 5.79 corresponding to i? = 0.05. We report our expression 
(I4.24P and the field-theoretical prediction of Refi^. For R = 0.03, our result is in very good agreement with 
the field-theoretical one. On the other hand, for R = 0.05, we observe significant differences for > 5. In 
any case, these discrepancies are within the 5% error we expect on our extrapolations. Indeed, for $p = 5 
we predict Zqs = 15.0 and Z{z = 5) = 12.9 to be compared with the field-theory results Zqs = 14.5 and 
Z{z = 5) = 12.7; for $p = 10 we obtain Zqs = 35.9 and Z{z = 5) = 28.5 to be compared with Zqs = 34.2 
and Z{z = 5) = 28.4. In all cases differences are less than 5%. Finally, note that both field theory and 
our results predict Z to be significantly lower than the good-solvent value Zqs for > 1 as soon as R is 
different from zero. 

A phenomenological expression for Z as a function of c and R — I — vf/vl/* is also reported in Refs.^§.i'2^. 
For ^ 1 there is good agreement'^^ with our results for all values of R, differences being less than 1%. 
On the other hand, significant differences are observed for larger values of $p.— For instance, for R = 0.03 
and ^p = 10 we predict Z/Zgs — 1 = —0.118, Ref<^ gives —0.113, in substantial agreement with our result, 
while the expression reported in Refi^ gives —0.060, which differs by a factor of 2. 

It is also interesting to compare the results for Rg/Rg. In Fig. [Tl] we report this ratio for $p = 1 as a 
function of z, together with the one-loop prediction of Ref.^ . In the good-solvent regime the difference is quite 
large. The poor behavior of the field-theoretical expressions in the dilute regime can be explained by looking 
at the virial expansion of Rg for z — + 00. In the good-solvent regime, the Monte Carlo simulations of RefJ^ 
give Slg « -0.315, S^g « -0.09, while Ref,^ predicts Sl g « -1.19, S^g « 5.74. One-loop perturbation 
theory renormalized as in Refi^ seems to be unable to reproduce correctly the polymer size as a function of 
$p, at variance with what occurs for Z. This is not totally surprising since the nonuniversal parameters that 
enter in the perturbative predictions (cq and no in the notation of Ref;^) were tuned to reproduce accurately 
the thermodynamic behaviori^ For instance, a different resummation of one-loop perturbation theory (see 
Sec. 6 of Refi^) gives Sl g w -0.145, 63 ,, ~ 0.0738, which are much closer to the Monte Carlo results. 

V. CONCLUSIONS 

In this paper we have determined the explicit TPM expressions for several quantities which characterize 
polymer solutions in the dilute regime. First, we have determined the universal constants An{z) for n = 
2, 3, 4. This allows us to obtain precise predictions for the osmotic pressure in terms of the polymer packing 
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fraction <i>p in the dilute regime $p < 1. Then, we have computed the swelhng factors for three different 
radii that characterize the polymer size. Finally, we have studied their concentration dependence. 

The expressions we have determined in this paper can be used in two different contexts. First, they 
provide expressions that may be used to fit the experimental data outside the universal (large- iV) regime. 
Note that the range of values of N in which a given system approximately behaves as predicted by the TPM 
expressions is nonuniversal, and thus in some cases the agreement is only at the level of the leading behavior, 
while in some others it may cover a significant range of polymer lengths. For instance, in some systems 
(e.g., in PMMA in chloroform or nitroethane at 20 °C) the interpenetration ratio approaches the universal 
value from above, i.e. > for large N (in the terminology of Ref.- , these systems are "strong-coupling 
systems" ) . This type of behavior cannot be described by the TPM expressions we have derived here which 
predict that the approach is always from below. In some other systems ( "weak- coupling systems") instead, 
the TPM expressions describe quite well the experimental behavior (for instance, polystyrene in cyclohexane 
or transdecaline, see chapter 15 in Rcf.^). Note also that 5" is generically a decreasing function of the 
temperature (at T — > Tg, \1/ — > 0) and thus a strong-coupling system becomes a weak-coupling system as T 
is lowered and thus eventually one can use the TPM to interpret the experimental behavior. 

The use of the TPM to describe the nonasymptotic behavior for finite values of N is mainly phenomeno- 
logical and rigorously justified only in systems in which the persistence length is larger than the typical 
monomer sizeJ^ A second use of the TPM expressions is in the description of the crossover to the 9 point. 
As we have explained in the introduction, the crossover functions defined in Eq. p.l[l can be computed in 
the TPM, by identifying z with {T — Te)N^^'^ (In N)~'^^^^ (modulo a normalization multiplicative constant). 
However, this identification is valid only close to the 9 point, since it assumes T — <C 1. To avoid this 
limitation, one can proceed as suggested in Refs3i2S, i.e., one can parametrize the crossover in terms of a 
physical variable. For instance, one can use the interpenetration ratio 4". Eq. (|l.ip can then be written as 

0{T,N,c) = aiOGiN,c)foi-^,cRl). (5.1) 

The quality of the solution is now characterized by ^E" that varies between (poor solvent) and (good 
solvent). In Fig. [T^lwe report the quantities computed above as a function of ^'/^E'*. Note that both A„ and 
Sn,# are small up to \E'/vE'* « 0.3 and also the swelling factors a do not change significantly in this range. 
This means that, for < 0.35** w 0.08, polymers behave approximately as Gaussian coils. Of course, as 
T — > three-body forces become increasingly important and thus tricritical corrections should be included. 
In the opposite range ^' > 0.08 tricritical effects can be neglected (see the numerical data in Ref.^) and one 
can use the TPM expressions to describe the polymer behavior. 



The authors thank Tom Kennedy for providing his efficient simulation code for lattice self-avoiding walks. 



APPENDIX A: PERTURBATIVE CALCULATIONS 



In this appendix we report the one-loop TPM expressions for the quantities reported in this paper. We 
use the general results of Ref."^-. For the osmotic pressure we start from the one-loop expression 

" ln[l + 2cur^ >{p)] + u / -— r , 



where u is the coupling constant and n is the polymer length. They are related to z and Rg by 

z = (27r)-3/2«ni/2 Rl^^ + 0{u). (A2) 

2 

The function T'^^\p) is the Debye function: 

4r(2)(p) = ^ - ^(1 - e-^^"/2). (A3) 



Expanding in powers of c we obtain A. 



4- 



1 31 072 

AAz) = -7777^(14075 + 12624V2 - lM&sVi)TT^''^ z^ + O(z^) = -29883. Iz'' -f Oiz^). (A4) 
45045 
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We give also numerical values for the leading TPM contributions to the following virial coefficients: 



(A5) 
(A6) 
(A7) 
(A8) 
(A9) 
(AlO) 

We can also compute the radius of convergence of the virial expansion. A simple analysis of the integral (|Aip 
shows that the singularity in the complex c-plane that is closest to the origin corresponds to cun^ = — 1. This 
implies An/An-i = —Sn^^^z asymptotically, and that the virial expansion converges for |<i>p| < l/(87r'^/^z) 
in the limit z 0. 

Using the expressions reported in Refi^ we obtain for the radius of gyration: 



Asiz) 


= 932283z^ 


Ae{z) 


= -3.13006- 10^ z^ 


Aviz) 


= 1.10136- 10^ 


Asiz) 


= -4.00539- 10l°z^ 


Agiz) 


= 1.49307- 10^2 


Aioiz) 


= -5.67361 - 10^^z^° 



Sl,g(z 

S2.,g{z 
Ss,g{z 



64 

3465 
2048 

405405 



(1365 - 1028\/2)7r3/2^2 ^ „9. 13421^2^ 

(85013- 115408\/2 + 45684V3)7r3z3 = 145.4282^ 



= -3113.85z^ 



Si,g{z) 

Sb,g{z) 

For the end-to-end distance we obtain analogously: 



78503.4z^ 
-2.19663- 10^ z^. 



Sl^e{z 

S2AZ 
Sz,e{z 



— (103- 76V2)7r3/2z2 
315^ ' 
4096 



40.1374z^ 



(1427 - 3248V2 + IS'i&^M)!^^ z'^ = 167.1322^ 



SsAz) 

Finally, for the hydrodynamic radius we use the representation 



10395 
-3654.20z^ 

93359. lz^ 

-2.6358- lO*' 



= 47r 



d^'q F{q) 
(27r)3 g2 ■ 



where F{q) is the form factor 



■(r-r,) 



(All) 

(A12) 

(A13) 
(A14) 
(A15) 

(A16) 

(A17) 

(A18) 
(A19) 
(A20) 

(A21) 
(A22) 



which is normalized so that F{q) = 1 — + 0{q'^). For the density corrections, we obtain at leading order 



Sim{z) = 3.32047z2 -I- 0(z3) 



S'2,//(z) 

Finally, we report the swelling factors: 



-44.8425z-^ + 0(z^ 



(A23) 
(A24) 



an 



1-f^ 
V 16 



4- ^log ^^)z + 0(z2) 



(A25) 
(A26) 
(A27) 
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Additional terms for are reported in Refi 
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TABLE I: Estiiriatos of A2(z) for five (iifforout values of z and different Nn^in- 

Z A'l.nii = 100 A'mm = 2:)() Aunu = ■")()() A'mjn = 1000 

Zl = 0.056215 0.99241(41) 0.99212(59) 0.99257(98) 0.9906(19) 

Z2 = 0.148726 1.97964(79) 1.9796(12) 1.9782(18) 1.9831(35) 

Z3 = 0.321650 2.9646(11) 2.9649(16) 2.9621(27) 2.9642(48) 

Z4 = 0.728877 3.9469(14) 3.9452(21) 3.9433(34) 3.9498(65) 

Z5 = 2.508280 4.9264(15) 4.9221(23) 4.9147(36) 4.9169(65) 



TABLE II: TPM estimates of several quantities for five different values of z. zi, Z2, Z3, za, and as are reported in the 
text. 





Zl 


Z2 


Z-A 


Za 




^2 


0.99257(98) 


1.9782(18) 


2.9621(27) 


3.9433(34) 


4.9147(36) 




0.08486(78) 


0.6061(30) 


1.8435(82) 


4.021(13) 


7.243(22) 


Aa 


-0.095(11) 


-0.857(78) 


-2.77(30) 


-5.51(79) 


-10.3(1.8) 


Age 


0.166145(65) 


0.165517(62) 


0.164591(65) 


0.163272(62) 0.161514(55) 


AgH 


1.50553(84) 


1.50964(90) 


1.51731(78) 


1.52597(74) 


1.54165(69) 


Sl,g 


-0.01711(66) 


-0.0630(11) 


-0.1278(18) 


-0.1970(24) 


-0.2747(25) 


'S'2,9 


0.0073(12) 


0.0309(42) 


0.051(10) 


0.082(17) 


-0.002(25) 


^l.e 


-0.01931(71) 


-0.0718(13) 


-0.1469(20) 


-0.2340(27) 


-0.3301(30) 


S2,e 


0.0097(14) 


0.0381(53) 


0.083(12) 


0.124(21) 


0.035(33) 


Sl^H 


0.00527(36) 


0.02079(75) 


0.0405(11) 


0.0602(13) 


0.0771(15) 


S2,H 


-0.00197(84) 


-0.0086(33) 


-0.0041(66) 


0.026(12) 


0.028(17) 


al 


1.06872(81) 


1.16553(82) 


1.31167(92) 


1.5683(12) 


2.2283(15) 




1.06545(71) 


1.15751(78) 


1.29540(89) 


1.5362(11) 


2.1595(13) 




1.03127(44) 


1.07196(44) 


1.12878(42) 


1.22196(40) 


1.43419(43) 



TABLE III: TPM interpolation formulas. 

A2{z) 47:^/^2(1 + 19.1187Z + 126.783^^ + 331.99^^ + 268.96z*)-^''* 

Aaiz) 1100.7z^(l + 23.1258Z + 195.358z^ + 594.386z^ + 541.9062"')"^/'* 

Aa{z) -29883.1z''(l + 17.4354z + 135.8532^ + 437A09z'^)-^^'^ 

Ageiz) 1/6 - 0.05714292(1 + 1455.412 + 47738.22^ + 584.5952^ + 5025.992^')-^/* 

AgH{z) 1.50451 + 0.02882522(1 - 1.644362 + 1.577772^ + 0.05350932^)"^/^ 

Si,g{z) -9.13421(1 + 0.01820932)2V(1 + 10.89442 + 28.83132^ + 0.527692^) 

52,9(2) 145.428(1 - 0.398682)2-^(1 + 25.O8O62 + 92.44152^ + 131.1642=^)-*/^ 

Si,e{z) -10.1374(1 + 0.016409l2)2V(l + 10.81712 + 26.19772^ + 0.4317312^) 

5'2,e(2) 167.132(1 - 0.1618192)2=^(1 + 24.65082 + 52.93132^ + 178.0512^)-'*/^ 

51, Hiz) 3.320472^(1 + 23.1962 + 71.94492^ + 257.6792=^ )-2/=^ 

52, Hiz) -44.8425(1 - 2.75132)2^(1 + 28.53272 - 37.08182^ + 350.1012^)"''/^ 
ag{z) (1 + 10.9288^; + 35.18692^ + 30.4463a;^)°°^*^*®'' 

ae{z) (1 + 11.41812 + 39.76612^ + 42.82572^ )°°5*='8«'^ 

OHiz) (1 + IO.435I2 + 29.76932^ + 16.89092^)""^*^*®'' 
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FIG. 1: The invariant ratio A2 as a function of for the five values of z reported in the text. We also report 

the linear extrapolation as determined by the fit of all data. 




FIG. 2: The ratios 6R^/N and Rh / [kuVN) (they converge to Qg and an, respectively) as a function of A'^ for 
the five values of z reported in the text. We also report the linear extrapolation as determined by the fit of all data. 
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FIG. 3: The TPM functions for the virial coefficients. We also show the Monte Carlo results reported in Table Ull 
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FIG. 4: The TPM functions for the density coefficients Sn.g, S„,e, and Sn,H- We also report the Monte Carlo results 
reported in Table |II1 
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FIG. 5: The TPM functions for the swelling factors. Together with our predictions we also report the results of 
Ref.^ (SCH) and of Ref."2^ (BN) for ag{z). Squares, triangles, and circles correspond respectively to the Monte Carlo 
results for Ofe, Qg, and an (see Table . 




FIG. 6: The compressibility factor vs $p for several values of z in the dilute region. 




FIG. 7: The compressibility factor vs $p for several values of z in the semidilute region. 




FIG. 8: Ratio R^/Rg for the radius of gyration for several values of z in the semidilute region. 




FIG. 9: The compressibility factor Z at <l?p = 1 as a function of z. We report our results ("our") and the field- 
theoretical ones reported in Ref.» ("SCH"). 




FIG. 10: Plot of Z/Zqs — 1 versus $p. Zqs is the compressibility factor in the good-solvent regime, while 
corresponds to solutions with two different values of i? = 1 — ^'/$*. We report the result (|4.24|l ("our") and 
field-theoretical one reported in Ref.- ("SCH"). 
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FIG. 11: The ratio Rg/R^ at "I>p = 1 as a function of z. We report our results ("our") and the field-theoretical ones 
reported in Ref.^ ("SCH"). 
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